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Abstract: We present the CMB anisotropy induced by the non- linear perturbations in 
the massive neutrino density associated to the non-linear gravitational clustering. We 
show that the non-linear time varying potential induced by the gravitational clustering 
process generates metric perturbations that affect the time evolution of the density fluc- 
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£N) , tuations in all the components of the expanding Universe, leaving imprints on the CMB 

anisotropy power spectrum at subdegree angular scales. For a neutrino fraction in agree- GO 
incut with that indicated by the astroparticle and nuclear physics experiments and a 
cosmological accreting mass comparable with the mass of known clusters, we find that 
CMB anisotropy measurements with Planck angular resolution and sensitivity possi- 
bly combined to other precise cosmological observations will allow the detection of the 
q dynamical, linear and non-linear effects of the neutrino gravitational clustering. 

h' 

CO 

> 

s_i ■ 1. Introduction 
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The atmospheric neutrino results from the Super-Kamiokande and MACRO ||] exper- 
iments indicate that neutrinos oscillate, these data being consistent with <-» v T oscil- 
lations. The small value of the difference of the squared masses (5 x 10~ 4 eV 2 < Am 2 < 
6 x 10 _3 eV 2 ) and the strong mixing angle (sin 2 28 > 0.82) suggest that these neutrinos are 
nearly equal in mass as predicted by many models of particle physics beyond the standard 
model. Also, the LSND experiment || support <-> v e oscillations (Am 2 < 0.2eV 2 ) and 
other different types of solar neutrino experiments suggest that u e could oscillate to 
a sterile neutrino v e <-> v s (Am 2 ~ 10 _5 eV 2 ). The direct implication of neutrino oscilla- 
tions is the existence of a non-zero neutrino mass in the eV range, and consequently a not 
negligible hot dark matter (HDM) contribution £l u ^ to the total energy density of the 
Universe. 
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We study the cosmic microwave background (CMB) secondary anisotropics induced 
by the non-linear perturbations in the massive neutrino density associated to the non- 
linear gravitational clustering. The extent to which the massive neutrinos can cluster 
gravitationally depends on their mass and the parameters of the fiducial cosmological 
model describing the present Universe. 

The linear perturbation theory describes accurately the growth of density fluctuations 
from the early Universe until a redshift z ~ 100. The solution involves the integration 



of coupled and linearized Boltzmann, Einstein and fluid equations [28] that describes the 
time evolution of the metric perturbations in the perturbed density field and the time 
evolution of the density fields in the perturbed space-time for all the relevant species (e.g., 
photons, baryons, cold dark matter, massless and massive neutrinos). At lower redshifts the 
gravitational clustering becomes a non-linear process and the solution relies on numerical 
simulations. 

Through numerical simulations, we compute the CMB anisotropy in the non- linear stages 
of the evolution of the Universe when clusters and superclusters of galaxies start to form 
producing a non-linear gravitational potential varying with time. By using a standard 
particle-mesh method we analyze the imprint of the dynamics of the neutrino gravitational 
clustering on the CMB anisotropy power spectrum in a flat ACHDM model with neutrino 
fractions /„ = n u /(Q b + O cdm )=0.06, 0.11, 0.16 corresponding to n„=0.022 (m„=0.78eV), 
0.037 (m„=1.35 eV), 0.053 (m !/ =1.89eV), assuming three massive neutrino flavors. This 
model is consistent with the LSS data and the WMAP anisotropy latest measurements 
JH, § allowing in the same time a pattern of neutrino masses consistent with the results 
from neutrino oscillation and double beta decay experiments. 
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2. The neutrino gravitational infall 

In the expanding Universe, neutrinos decouple from the other species when the ratio of their 
interaction rate to the expansion rate falls below unity. For neutrinos with masses in the eV 
range the decoupling temperature is Tp ~ lMeV, occurring at a redsfit Z£> ~10 10 0. At 
this time neutrinos behave like relativistic particles with a pure Fermi-Dirac phase-space 
distribution: 

f v (q,a) = E J Ti/+1 , K = vV + a?ml , (2.1) 

where q is the neutrino comoving momentum, q = ap, p being the neutrino 3- vector mo- 
mentum, E v is the energy of neutrino with mass m u and a = 1/(1 + z) is the cosmic scale 
factor evolving with the time, t (clq=1 today). 

As neutrinos are collisionless particles, they can significantly interact with photons, baryons 
and cold dark matter particles only via gravity. The neutrino phase space density is con- 
strained by the Tremaine & Gunn criterion l| that put limits on the neutrino energy 
density inside the gravitationally bounded objects: in the cosmological models involving 
a HDM component (the CHDM models) the compression fraction of neutrinos through a 
cluster /(r) = p u /Pcdm (where r is the cluster radius) never exceeds the background ratio 
^W^edm ||- Because the formation of galaxies and clusters is a dynamical time process, 
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Figure 1: The evolution with the redshift of the projected mass distributions of cold dark matter 
plus baryons (upper row) and neutrinos (lower row) obtained from numerical simulation of 128 3 
cold dark matter plus baryons (the total mass of 8 x 1O 16 M ) and 10 x 128 3 neutrinos (the total 
mass of 4.8 x 10 18 M Q ) in a box of size 128 Mpc, for the ACHDM model with the neutrino fraction 
f u = 0.06 (J2i m vi ~ 0.7 eV). [Units of axes are in Mpc]. 



the differences introduced in the gravitational potential due to neutrino gravitational clus- 
tering generate metric perturbations that affect the evolution of the density fluctuations 
of all the components of the expanding Universe. Fig. 1 presents the evolution of the pro- 
jected mass distributions of cold dark matter plus baryons and neutrinos obtained from 
numerical simulations at few redshift values z. One can see that neutrinos are accreted 
by the cold dark matter and baryons, contributing in dynamic way to the gravitational 
clustering process. 

Neutrinos cannot cluster via gravitational instability on distances below the free-streaming 
distance Rf s 11, |l2|| . The neutrino free-streaming distance is related to the causal co- 



moving horizon distance t](a) through [p~3|: 



Rf s (a) 



7](a) 
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a 2 H( a y 



(2.2) 



where a nr is the value of the scale factor when massive neutrinos start to become non- 
relativistic (a nr = (1 + z nr )~ x « 2>kBT v $/m y c 2 ) and H(a) is the Hubble expansion rate: 
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Figure 2: Panel a): dependence of the neutrino free-streaming distance Rf s on the scale factor. 
We show also a specific scale, A=20/i _1 Mpc, constant in comoving coordinates (horizontal dashed 
line) and the evolution with the scale factor of the causal horizon distance for our cosmological 
models. Panel b): dependence on the scale factor of the mass M(Rf s ) of the perturbation at the 
scale Rf s . We show also the mass of the perturbation at the scale A (horizontal dashed line) and 
indicate the typical mass ranges for galaxies (GAL), groups (GR), clusters (CL) and superclusters 
(SCL). 



In the above equation G is the gravitational constant, D, m = f2& + SI dm + Sl u is the matter 
energy density parameter, 0^, fi C( i m , Q v being the energy density parameters of baryons, 
cold dark matter and massive neutrinos, £l r is the radiation energy density parameter that 
includes the contribution from photons and relativistic neutrinos , is the vacuum (or 
cosmological constant) energy density parameter, Qk = 1 — U m — CIa is the energy density 
parameter related to the curvature of the Universe. 

Rf s defines the minimum linear dimension that a neutrino perturbation should have in 
order to survive the free-streaming. In the spherical approximation, the minimum comoving 
mass of a perturbation that should contain clusterized neutrinos, corresponds to Q] 

M(R fs ) = ^R%p m « 1.5 x 10 11 (O m / l 2 )(i? /s /Mpc) 3 / l - 1 M , 

where Q m is the matter energy density parameter. 

We show in Fig. 2 the dependence of the causal horizon distance 77(a), the neutrino free- 
streaming distance Rf s , [panel a)] and of the mass M(Rf s ) [panel b)] on the cosmic scale 
factor. The cosmological model is the ACHDM model with different neutrino fractions 
f v . One can see that at early times, when neutrinos are relativistic, the free-streaming dis- 
tance is approximately the causal horizon distance. After neutrinos become non-relativistic 
(a nr ~ 10~ 4 for our cosmological models) the free-streaming distance decreases with time, 
becoming smaller than the causal horizon distance. The time behaviors of Rf s and M{Rf s ) 
show that neutrino can cluster gravitationally on increasingly smaller scales at latter times. 
If the causal horizon 77(a) is large enough to encompass the wavelength A, the neutrino grav- 
itational infall perturbs the growth of the perturbations for this mode, leaving imprints in 
the CMB angular power spectrum. Perturbations on scales A < Rj s (k > kf s ) are damped 
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due to the neutrino free-streaming while the perturbations on scales A > Rf s (k < kf s ) are 
affected only by gravity. In the Newtonian limit, the neutrino gravitational clustering can 
be described as a deviation from the background by a potential $ given by the Poisson 
equation: 

V 2 $(r,a) = 4irGa 2 Pm (a)5 m (r,a) , (2.4) 

where f is the position 3-vector, p m (a) is the matter density and 5 m (f,a) is the matter 
density fluctuation; 5 m = (pbdb + Pcdm^cdm + Pv^v)/Pm, where p c d m , Pb, Pv and 5 c d m , 4, 
5 U are the density and density fluctuations of cold dark matter particles, baryons and 
neutrinos. 

The equations governing the motion of each particle species (cold dark matter plus baryons 
and neutrinos) in the expanding Universe are given by pi 



^ = -aH(a)VQ, % = q '(a 3 H (a))- 1 , (2.5) 
da da 

where q is the comoving momentum and H(a) is given by the equation (2.3). 
The Newtonian description given by the equations (2.4) and (2.5) applies in the limit of 
the week gravitational field if, at each time step, the size of the non- linear structures is 
much smaller than the causal horizon size (the background curvature is negligible). 

The cosmological models involving massive neutrinos show a characteristic scale de- 
pendence of the perturbation growth rates [pi], 17, 18 1. 



We evolve the system of baryons plus cold dark matter particles and neutrinos accord- 
ing to the equation (2.5) for the non-linear scales involved in the computation of the CMB 
anisotropy (0.06Mpc -1 < k <0.52Mpc _1 ), starting from the beginning of the non-linear 
regime of cold dark matter plus baryons component. 

The initial positions and velocities of neutrinos and baryons plus cold dark matter particles 
can be generated at each spatial wave number k from the corresponding matter density 
fluctuations power spectra at the present time by using the Zel'dovich approximation ]l9| . 
The matter power spectra was normalized on the basis of the analysis of the local cluster 
X-ray temperature function [2C[. We performed simulations with 128 3 cold dark matter 



plus baryon particles and 10 x 128 3 neutrinos. The neutrinos and the baryons plus cold 
dark matter particles was randomly placed on 128 3 grids, with comoving spacing vq of 
0.5 h _1 Mpc. The high number of neutrinos and this comoving spacing ensure a precision 
high enough for a correct sampling of the neutrino phase space distribution [18]. 



According to the Zel'dovich approximation, the perturbed comoving position of each par- 
ticle r(r*o, a) and its peculiar velocity v(ro, a) are related to the fluctuations of the density 
field 5p(ro, o, k) through: 

f(f ,k,a) = f + D(k,a)d(f ) , v(f ,k,a) = D(k,a)d(f ) , (2.6) 
Vd(r ) = D~ 1 (k,a)5p(r , k,a) , 

where ro is the coordinate corresponding to the unperturbed comoving position, d(fo) is 
the displacement field and D{k, a) is the growth function of perturbations corresponding 
to each cosmological model. 
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Figure 3: The scale dependence of z n i on fc for neutrinos (panel a) and baryons plus cold dark 
matter particles (panel b). 



At each wave number k used in the computation of the CMB anisotropy we compute the 
perturbed particles comoving positions and peculiar velocities at the beginning of the non- 
linear regime a n i, by using the set of equations (2.4)-(2.5). We assign to each particle a 
momentum according to the growth function, when the power of each mode is randomly 
selected from a Gaussian distribution with the mean accordingly to the corresponding 



power spectrum |2l| , 22, \Z3^ . In the computation of the set of equations (2.6) we consider 
only the growing modes, the non-linear power spectra up to k max = 6.28 h Mpc -1 , and 
neglect the contribution of the redshift distortions. 

We show in Fig. 3 the dependence on the spatial wave number k of the redshift z n i = 
^/ a ni ~ 1 f° r each component. One can see from Fig. 3 that neutrinos (panel a) enter in the 
non- linear regime later than cold dark matter particles and baryons (panel b). Thus, the 
neutrino halo of the cluster starts to form after the cold dark matter plus baryon halo is 
advanced in the non-linear stage, causing the accretion of neutrinos from the background. 
At each spatial wavenumber k we evolve the particles positions and velocities according 
to the set of equations (2.4)-(2.5). We start this process from the scale factor a^j m at 
which cold dark matter particles plus baryons start to enter in the non-linear regime. At 
each time step, the density on the mesh is obtained from the particle positions using the 
Cloud-in-Cell method and equations (2.6) are solved by using 7-point discrete analog of the 
Laplacian operator and the FFT technique. The particle positions and velocities are then 
advanced in time with a time step da required by the computation of the CMB anisotropy 
power spectra. The system of particles was evolved until the scale factor a s t when it reaches 
its virial [p4| equilibrium. 



3. Imprints of neutrino gravitational clustering at Planck angular scales 

As the anisotropy produced by the non-linear density perturbations depends on the time 
variations of the spatial gradients of the gravitational potential produced by different com- 
ponents (cold dark matter, baryons, neutrinos), we calculate the CMB anisotropy in the 
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Figure 4: Panel al): time evolution of the energy density perturbations of the different components 
as computed by including linear and non-linear effects of neutrino gravitational clustering (solid 
lines) and by neglecting the non- linear aspects of neutrino gravitational clustering (dashed lines): 
cold dark matter (cdm) , baryons (bar) , massive neutrinos (y) , and massless neutrinos plus photons 
(rad). Panel a2): the same as in panel a), but for the time evolution of the gravitational field 
[k^.OGMpc- 1 and /„ = 0.06]. 



presence of the gravitational clustering by using N-body simulation in large boxes with the 
side of 128 Mpc, that include all non-linear scales used in the computation of the CMB 
anysotropy power spectrum from \ m in ~ 12Mpc (k max « 0.52Mpc _1 ) to X ma x ~ 110 Mpc 
(k m i n ~ 0.06Mpc _1 ), taking into account the time evolution of all non-linear density per- 
turbations influencing the CMB power spectrum (see also Fig. 1). One should note that 
Xmax corresponds to the comoving horizon size at the matter-radiation equality for our 
cosmological models (X eq ~ 16 Q~ /i —2 Mpc). The non-linear structures are assumed to be 
formed by two components: cold dark matter plus baryons and neutrinos in the form of 
three massive neutrino flavors, both components evolving in the gravitational field created 
by themself. For the purpose of this work we neglect the hydrodinamical effects [25]. The 
neutrino gravitational clustering can affect both the homogeneous and the inhomogeneous 
components of the gravitational field. The changes in the homogeneous component of the 
gravitational field are determined by the changes of the energy density of neutrinos and 
cold dark matter particles plus baryons. They affect the Hubble expansion rate, the sound 
horizon distance and the neutrino free-streaming distance. The changes in the inhomoge- 
neous component of the gravitational field are determined by the changes in the energy 
density for all matter components and the changes in the neutrino phase space distribution 
function. They affect the growth of the energy density perturbations of cold dark matter, 
baryons, photons, massive and massless neutrinos. Panel al) of Fig. 4 presents the evolu- 
tion with the scale factor of the energy density perturbations of different components in 
the non-linear regime, for the mode k = 0.06 Mpc -1 (solid lines). For comparison, we plot 
also (dashed lines) the energy density perturbations of the different components obtained 
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Figure 5: The imprint of the neutrino linear and non-linear gravitational clustering on the CMB 
anisotropy power spectrum expressed in terms of AT/T obtained for the filtering perturbation with 
the mass M. We report in the panels the richness of these perturbations from the bottom to the top 
according to the increasing of the power at multipoles about I ~ 750. In each panel, the dashed line 
corresponds to the fiducial ACHDM cosmological model, without including the non-linear effects of 
neutrino gravitational clustering. 

for the same mode k by neglecting non-linear aspects of neutrino gravitational clustering 
(linear regime). Panel a2) of Fig. 4 presents the evolution with the scale factor of the scalar 
potential of the conformal Newtonian gauge line element, that plays the role of the grav- 



itational potential in the Newtonian limit [26, 27], by including (solid line) or not (dashed 
line) the non-linear effects of neutrino gravitational clustering (for the transformation re- 
lation between the scalar potentials of the synchronous gauge and conformal Newtonian 
gauge see equation (18) from |28|| ). As we have shown before, the difference in the evolution 
of a perturbation mode k depends on how this mode relates to the neutrino free-streaming 
wave number kf s . Considering that our simulation at each time step is a sample of the 
evolution of the matter in the non-linear regime, we study the imprint of the gravitational 
clustering on the CMB anisotropy power spectrum by smoothing the density field obtained 
from simulation at each time step with a filter with the scale Rf s corresponding to the 
cluster mass value M(Rf s ). For each non- linear mode k only the perturbations with the 
mass M < M(Rf s ) are taken into account for the computation of the CMB anisotropy 
power spectrum. Fig. 5 presents some of our computed CMB anisotropy power spectra 
obtained when different filtering mass values M(Rf s ) are considered. It is usual to use the 
Coma cluster as the mass normalization point (Mcoma = 1-45 x 10 15 h _1 MQ); for the Coma 
cluster we assume a richness *4coma=106. According to |P9| i the relation between the mass 
of the perturbation and the richness A of the corresponding cluster can be written in the 
form 



M = M Co 



A 



-4coma 



1.45 x 10 



15 



A_ 
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hr x M, 







By comparing the angular power spectra obtained including or not the non-linear effects of 
neutrino gravitational clustering, we find a decrease of the CMB angular power spectrum 
induced by the neutrino non-linear gravitational clustering of AT/T ~ 10 -6 for angular 
resolutions between ~ 4 and 20 arcminutes, depending on the cluster mass and neutrino 
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fraction f u . 

Clearly, new high sensitivity and resolution anisotropy experiments will have the capa- 
bility to detect the neutrino gravitational clustering effect. In particular, the instruments 
on-board the ESA Planck satellite 1 will measure the CMB angular power spectrum with 
very high sensitivity up to multipoles i ~ 1000 — 2000 with a stringent control of the sys- 
tematic effects. Fig. 6 compares Planck and WM AP 2 performances. The CMB angular 
power spectrum is reported without beam smoothing and by taking into account the beam 
window functions of several Planck frequency channels and of the highest WMAP fre- 
quency channel (which is very close to that of the LFI 70 GHz channel). The corresponding 
angular power spectra of the residual nominal white noise (i.e. after the subtraction of the 
expection of its angular power spectrum) are also displayed. Of course, binning the power 
spectrum on a suitable range of multipoles, as usual at high £, will allow to recover the 
CMB power spectrum also at multipoles higher that those corresponding to the crossings 
between the noise and CMB power spectra reported in the figure 3 . 

The characteristic angular scale left by the the neutrino gravitational clustering on the 
CMB anisotropy power spectrum is given by 

Q = R fs 

% - 77(a) ' 

where Rf s is the scale of the filtering perturbation with the mass M(Rf s ), 77(a) is the 
particle horizon distance at the time at which the non-linear perturbation mode k cross 
the horizon and 770 is the particle horizon at the present time. Fig. 7 (left panel) presents 
the evolution of the characteristic scale 6 and of the corresponding multipole order of the 
CMB anisotropy power spectrum with the mass M(Rf s ). 

Fig. 7 (right panel) presents few confidence regions of the f u — M parameter space 
that can be potentially detected by the Planck surveyor by using the CMB anisotropy 
measurements in the presence of the gravitational clustering, under the hypothesis that 
the other cosmological parameter can be measured with other observations (grey regions) 
and by using the Planck data to jointly determine f u , M and the other cosmological 
parameters (solid lines). We consider for this computation only the Planck "cosmo- 
logical" channel between 70 and 217 GHz, a sky coverage f s ^ y = 0.8 and neglect for 
simplicity the foreground contamination |fTq| . By assuming known the other main cos- 
mological parameters, we obtain a neutrino fraction /„ ~ 0.011 ± 0.007 for an accreting 
mass M ~ (8.2 ± 3.1) x 1O 14 /i _1 M0 (errors at 68% confidence level). By assuming known 
the other main cosmological parameters, we obtain a neutrino fraction /„ ~ 0.011 ± 0.007 
for an accreting mass M (8.2 ± 3.1) x lO 14 n _1 M (errors at 68% confidence level). 
PLANCK surveyor will have in principle the capability to measure the non-linear imprints 
of the neutrino gravitational clustering on the CMB anisotropy power spectrum for a neu- 
trino mass range in agreement with that indicated by the astroparticle and nuclear physics 

1 http: / / astro.estec.esa.nl/Planck 
2 http:/ /lambda. gsfc.nasa.gov 

3 At £ higher than ~ 1500 — 2000 the confusion noise from extragalactic source fluctuations dominates 
over the instrumental noise, while Galactic foregrounds are relevant at multipoles less than few hundreds. 
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Figure 6: Comparison between Planck and WMAP resolution and sensitivity. For each con- 
sidered frequency channel, the crossing between the CMB convolved angular power spectrum and 
the unsubtractable instrumental white noise angular power spectrum indicates the multipole value 
where the signal to noise ratio (£ by t) is close to unity. 

experiments and a cosmological accreting mass comparable with the mass of the known 
clusters. Of course, even with the high sensitivity and resolution of Planck it is hard 
to firmly constrain f u and M by jointly recovering the other cosmological parameters, a 
goal that can be achieved in combination with other precise cosmological information, such 
as galaxy large scale structure surveys, measures of element abundances from big-bang 
nucleosynthesis and Type la supernovae observations. 

4. Conclusions 

We study the CMB anisotropy induced by the non-linear perturbations in the massive 
neutrino density associated to the non-linear gravitational clustering. Through numerical 
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Figure 7: Left panel: evolution of the characteristic angular scale of the neutrino gravitational 
clustering with M(Rf s ) for different neutrino fractions. We report also the corresponding multi- 
pole orders I ~ B^ 1 of the CMB anisotropy power spectrum. Right panel: confidence regions of 
the /„ — M parameter space as potentially detectable by Planck by using the CMB anisotropy 
measurements in the presence of the gravitational clustering assuming known the other main cos- 
mological parameters (grey regions) or by jointly recover /„, M and the other main cosmological 
parameters (solid lines). 



simulations, we compute the CMB anisotropy angular power spectrum in the non-linear 
stages of the evolution of the Universe when clusters and superclusters start to form, pro- 
ducing a non-linear time varying gravitational potential. 

We found that the non-linear time varying potential induced by the gravitational clustering 
process generates metric perturbations that affect the time evolution of the density fluc- 
tuations in all the components of the expanding Universe, leaving imprints on the CMB 
anisotropy power spectrum at subdegree angular scales. The magnitude of the induced 
anisotropy and the characteristic angular scale depends on how each non-linear mode k of 
the perturbations relates to the neutrino free-streaming wavenumber kf s at each evolution 
time step. By smoothing the density field obtained from simulations with a filter with 
the scale corresponding to the cluster scale, we find an imprint on the CMB anisotropy 
power spectrum of amplitude AT/T ps 1CT 6 for angular resolutions between ~ 4 and 20 
arcminutes, depending on the cluster mass and neutrino fraction f u . 

This result suggests that the CMB anisotropy experiments with such levels of sensitivities 
and angular resolutions should detect the dynamical effect of the non-linear gravitational 
clustering. For a neutrino fraction in agreement with that indicated by the astroparticle 
and nuclear physics experiments and a cosmological accreting mass comparable with the 
mass of known the clusters, we find that CMB anisotropy measurements with Planck 
angular resolution and sensitivity in combination with other precise cosmological observa- 
tions will allow the detection of the dynamical, linear and non-linear, effects of the neutrino 
gravitational clustering. 
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